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^ . ABSTRACT 

m ■ 

Context. Despite the low cosmic abundance of deuterium (D/H ~ lO"""), large degrees of deuterium fractionation in molecules are observed in 
star forming regions with enhancements that can reach 13 orders of magnitude, which current models have difficulties to account for. 
Aims. Multi-isotopologue observations are a very powerful constraint for chemical models. The aim of our observations is to understand the 
processes forming the observed large abundances of methanol and formaldehyde in low-mass protostellar envelopes (gas-phase processes ? 
O I chemistry on the grain surfaces ?) and better constrain the chemical models. 

Methods. Using the IRAM 30m single-dish telescope, we observed deuterated formaldehyde (HDCO and D2CO) and methanol (CH2DOH, 
CH3OD, and CHD2OH) towards a sample of seven low-mass class protostars. Using population diagrams, we then derive the fractionation 
ratios of these species (abundance ratio between the deuterated molecule and its main isotopologue) and compare them to the predictions of 
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grain chemistry models. 

• • Results. These protostars show a similar level of deuteration as in IRAS 16293-2422, where doubly-deuterated methanol - and even 
, triply-deuterated methanol - were first detected. Our observations point to the formation of methanol on the grain surfaces, while formaldehyde 
, formation cannot be fully pined down. While none of the scenarii can be excluded (gas-phase or grain chemistry formation), they both seem to 
» I require abstraction reactions to reproduce the observed fractionations. 

a 
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1. Introduction 

In the last few years, observations of low-mass protostars have 
revealed unexpected high abundances of deuterated molecules, 
and particularly doubly-deuterated molecules. The discov- 
ery of an extremely large amount (D2CO/H2CO ~ 10%) 
of doubly-deuterated formaldehyde in the l ow-mass proto- 
star IR AS 16293 -2422 (hereafter IRAS 1 6293 JCeccarelli et all 
Il998 ») was followed by observations of this same molecule to- 
wards a large sample of low-mass protostars (Loinard e t all 
I2OO2I) . The fractionation ratios appeared to be similarly large 
towards all targeted low-mass protostars. The suggested inter- 
pretation was that such large deuteration is obtained during the 
cold and dense precollapse phase of the low-mass protostars: 
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highly-deuterated ices are very likely formed via active grain 
surface chemistry (Tielens 1983), stored in the grain mantles 
and eventually released in the gas phase during the collapse, 
when the heating of t he new ly-formed protostar evaporates the 
ices (Ceccarelli et al1 l200lh . A strong support to this scheme 
comes fro m the large abundanc e of D2CO observed in prestel- 
lar cores (jBacmann et al. hoof), and from the discovery o f the 
large fractio nation of in the same objects (Caselli e t alJ 
l2003t IVasteletalJ 120041) . While these last observations un- 
doubtly supported the idea that deuteration sets on just before 
the collapse, the leading chemical process was still largely un- 
proven. In this context, methanol is a key molecule. Indeed, it 
is believed to be a grain surface product as gas-phase processes 
are not efficient enough to acco unt for the large abundances 
observed in star-forming regions ( iHerbsi200 5). Methanol may 
be the last step of CO hydrogenation on the grain surfaces, 
after the formaldehyde formation. If formaldehyde formation 
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is also dominated by grain chemistry, the deuteration of both 
molecules is expected to be tightly linked. 



Three years ago, doubly-deuterated methanol was 
detected towards the low- mass protostar IRAS 16293 
(CHD2OH/CH3OH ~ 20%, I Parise et alJ l2002h . basically 
confirming the grain chemistry scheme. However, the 
CH^DOH/ CH^OD ratio was f ound to be unexpectedly large 
(20 + 14, IParise et all 12002') compared to the value of 3 
predicted by the models (Charnlev et al. 1997). The proposed 
hypothesis of a rapid conversion of CH3OD into CH3OH in the 
gas phase due to protonation reactions that would affect only 
the species for which the deuterium is bound to the very elec- 
tronegative oxygen ( Charnlev et alJll9 97) was also suggested 
by the model of Osa mura et al.l ll2004l) . The high observed 
fractionation of CH2DOH and CHD2OH was consistent with 
formation of methanol on the grain surfaces, but required an 
atomic D/H ratio in the gas-phase as high as 0.1-0.2 during 
the mantle formation. This v alue was cha l lengin g gas-phase 
models at that time, and led IParise et alJ ( l2002h to suggest 
that "a key parameter was missing" in the chemical schemes. 
Indeed, this key parameter was soon discovered. In molecular 
clouds, the main reservoir of deuterium is molecular HD. 
Some deuterium can be transferred from this reservoir to other 
molecules by the intermediate of the ion H2D^ that forms 
according to the exothermic reaction: HD H- H^ ^ H2 H- H2D^. 
Collision of CO with H2D^, followed by recombination will 
then form atomic deuterium. Depletion of CO — one of the 
main destruction agents of H2D^ — drives the fractionation of 
this key intermediary in the gas phase deuterium fractionation 
schemes all the way to HDJ and even DJ . The inclusion into 
the gas-phase schemes of multiply-deuterated isotopomers of 
H2D^ as new intermediate molecules for deuterium transfer 
from the HD main reservoir allowed to pr edict the atomic D/H 
ratio required by methanol observations ( Roberts et al. 2003). 
Triply-deuterated methanol was later detected in IRAS 16293 
(Parise et al. 2004), and the CH3OH column density was better 
evaluated by analyzing '^CH30H transitions. The observed 
CD3OH/CH3OH fractionation ratio, that was found to be 
consistent with the CH2DOH and CHD2OH fractionations, 
allowed to confirm the grain chemistry schem e. While much 
progr ess has been made through laboratory (e.g. lNagaoka et all 
I2OO5I) and theoretical studies to understand the observed high 
deuterium fractionation in IRAS 16293, observations of a 
larger sample are required to determine whether this is a 
common phenomena in protostellar environments. 



In this paper, we report singly (HDCO) and dou- 
bly (D2CO) deuterated formaldehyde and singly-deuterated 
methanol CH2DOH and CH3OD as well as doubly-deuterated 
CHD2OH observations towards a sample of low-mass proto- 
stars. This article is organized as follows: observations are pre- 
sented in section 2, analysis of the data is presented in section 
3, the results are discussed in section 4 and conclusions are 
drawn in section 5. 



2. Observations and results 

2.1. Observations 

Using the IRAM 30-meter telescope (Pico Veleta, Spain), we 
observed the five deuterated species HDCO, D2CO, CH2DOH, 
CH3OD and CHD2OH towards the six low-mass protostars 
NGC1333-IRAS4A, -IRAS4B, -IRAS2,L1448N,L1448mm 
and L 11 57mm. We also present observations of deuterated 
formaldehyde (HDCO and D2CO) towards L1527. All these 
sources, already studied by Maret et al. (2004), are Class 
protostars, i.e. in the early phase of the gravitational collapse. 
Table^summarizes the main characteristics of these sources. 

The observations were performed in September and 
November 2002, September 2003 and March 2004. For 
methanol, four receivers were used simultaneously at 3, 3, 1.3 
and 1.3 mm with typical system temperatures of about 100, 
150, 400 and 400 K respectively. These receivers were con- 
nected to the VESPA autocorrelator, divided in 6 units. For 
formaldehyde, four receivers were used simultaneously at 3, 
2, 1.3 and 1 mm, with typical system temperatures of about 
110, 220, 300, and 400 K, connected to VESPA. The tele- 
scope beam width varies between 30" at 83 GHz and 9" at 
276 GHz. All observations were performed using the position 
switching mode. The offset positions are summarized in Table 
[2 The pointing accuracy was monitored regularly on strong ex- 
tragalactic continuum sources and found to be better than 3". 
Our spectra for deuterated methanol were obtained with inte- 
gration times ranging from 280 to 420 minutes depending on 
the source. 

All fluxes are indicated in units of main-beam temperature, 
and the error bars were calculated as the quadratic sum of the 
statistical noise (noise rms of the data) and the calibration un- 
certainty. In order to account for the atmospheric calibration as 
well as for uncertainties in the band rejection, and taking into 
account that most of our observations were performed at high 
elevation, we adopted the following calibration uncertainties: 
5% for lines observed at frequencies lower than 130 GHz, 10% 
for frequencies between 130 and 260 GHz, and 15% for higher 
frequencies. 

All upper limits are given at a 3cr level: 

J^Tnibdv < 3cr(l -t- a) ^6v.Av 

where a is the calibration uncertainty indicated above, <t the 
noise rms of the observations, 5v the spectral resolution and 
Av the assumed linewidth (1.5 kms ' for D2CO, 3 kms"' for 
deuterated methanol), based on observations of detected lines. 

Observed fluxes are listed in Tables|6lE]|8land|9l Examples 
of observed spectra are shown in Fig.l7ll8ll9l [T0lfm and[T2l 

2.2. Results 

We detected the two deuterated formaldehyde isotopes (HDCO 
and D2CO) towards all the sources of our sample. Regarding 
methanol, only the three sources NGC1333-IRAS4A, 
-IRAS4B and -IRAS2 ha ve good eno ugh detections. This is 
consistent with the study bv lMaret et ali (.2005.) that shows that 
these three sources are indeed the brightest ones for CH3OH 
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Table 1. Parameters of the sources targeted in our deuterated formaldehyde and methanol observations. IRAS 16293-2422 is 
added for comparison. The offset used for the OFF position of the position switching mode are indicated in the last column. 



Source 


a (2000) 


S (2000) 


Region 


Distance 


T " 


M 


Lsmm/Lbol 


Tboi 


^LSR 


offsets 










(parsec) 


(Lo) 


(Mo) 


(%) 


(K) 


(kms-1) 


(arcsec) 


IRAS4A 


03:29:10.3 


31:13:32 


Perseus 


220 


6 


2.3 


5 


34 


+7 


(-150,240) 


IRAS4B 


03:29:12.0 


31:13:09 


Perseus 


220 


6 


2.0 


3 


36 


+7 


(-180, 260) 


IRAS2 


03:28:55.4 


31:14:35 


Perseus 


220 


16 


1.7 


< 1 


50 


+7 


(70, 180) 


L1448N 


03:25:36.3 


30:45:15 


Perseus 


220 


6 


3.5 


3 


55 


+5 


(300, 300) 


LI 448mm 


03:25:38.8 


30:44:05 


Perseus 


220 


5 


0.9 


2 


60 


+5 


(300, 300) 


LI 157mm 


20:39:06.2 


68:02:22 


Isolated 


325 


11 


1.6 


5 


60 


+5 


(150, 0) 


LI 527 


04:39:53.9 


26:03:10 


Taurus 


140 


2 


0.9 


0.7 


60 


+5 


(240, 0) 


IRAS 16293 


16:32:22.6 


-24:28:33.0 


p-Ophiucus 


160 


27 


54 


2 


43 


+4 





From'Andreetal.' ("200(f) and' CemisI il99ol) 
* From .lorgensen et aLl200Z) 



emission. CH2DOH was detected in all sources in which it was 
searched for, but only the low-lying transition was detected 
in the case of L1448N, L1448mm and L1157mm. CH3OD 
was detected only towards IRAS4A, IRAS4B, L1448mm and 
LI 157mm, with only one transition detected in the two last 
sources. Finally, doubly-deuterated methanol was detected 
towards IRAS4A, IRAS4B and IRAS2. Upper limits were de- 
rived for the other sources. Despite the substantial integration 
time (420 minutes for L1448N, 108 minutes for L1448mm and 
264 minutes for L 11 57mm), upper limits on the fractionation 
are not very significant, m ostly because of t he low CH3OH 
abundance in those sources dMaret et alJ2005h . 

The observed lines are relatively narrow, around 1 to 1.5 
kms ' for deuterated formaldehyde and up to 3 kms ' for 
deuterated methanol. 

In particular for formaldehyde, the linewidths for the 
deuterated species are smaller than for the main isotopomer 
(Maret et al. 2004) . Indeed, the low energy H2CO lines are 
pr obably contaminate d by an outflow contribution, as discussed 
in Maret eta l] |l2004 . The narrow lines emitted by the deuter- 
ated isotopomers of H2CO suggest that the emission is domi- 
nated in this case by the cold outer envelope of the protostar, 
as expected from such lines with relatively low upper ener- 
gies. This is in agreement with the observations of extended 
emission of D7CO in th e low-mass protostar IRAS 1 6293-2422 
dCeccarelli et al.l200ll) . This observation was interpreted as an 
indication that D2CO comes from the evaporation of CO-rich 
ices that evaporate around 20 K, i.e. at a lower temperature than 
polar ices. 

For methanol, the signal to noise ratio is too low to draw 
any firm conclusion, but the deuterated lines also seem to be 
narrower than the main isotopomer, suggesting that the emis- 
sion is also in this case dominated by the envelope contribution. 
However, the larger linewidths compared to formaldehyde as 
well as the higher rotational temperature (see next paragraph) 
suggests that it probes more deeply the envelope. 



3. Analysis 

The analysis based on an accurate model of infalling envelope 
of H2CO and CH3OH transitions towards the sample of low- 
mass protostars showed that the abundance o f those two species 
jumps in the inner warm part of the envelope ( Maret et al.'2 004 
2005). This abundance jump has been attributed to the evapora- 
tion of H2CO and CH3OH from polar ices, in the region where 
the temperature is higher than lOOK. In principle, a multifre- 
quency analysis could be done for the CH2DOH molecule, for 
several lines have been observed in this molecule. However, the 
analysis could only be done in the LTE approximation, for the 
collisional coefficients are not known. Besides, since we cannot 
perform the analysis for the other molecules, we would not be 
able to make any comparison anyway. We therefore decided not 
to implement this analysis in this article. We will, on the con- 
trary, use the rotational diagrams technique, which gives the 
column density averaged on the source extent for optically thin 
and LTE lines. Both conditions are likely correct in our case, 
first because we do not expect particularly large column den- 
sities, second because the critical densities' for these lines are 
around 10^ to 10^ cm""* which are about the densities around 
these sources jMaretetalJ2004l> . 

The lines have been observed at several different frequen- 
cies, hence with different spatial resolution. Indeed, the beam 
size of the IRAM 30 m telescope is 30" at 83 GHz and 9" at 
276 GHz. If the source is smaller than the beam size of the ob- 
servation, the derived column density must be corrected for the 
beam dilution. Unfortunately, we do not know for sure the size 
of the deuterated formaldehyde and methanol emission in the 
targetted protostars. In order to check whether and how much 
the uncertainty on the source size affects the derived column 
densities and fractionations, we discuss in the following para- 
graph a detailed study of the CH2DOH data on IRAS4A. 



' For example, for CH2DOH, the Einstein coefficients for the ob- 
served transitions lie in the 10"' to 10"^ s"' range. The collisional de- 
excitation rates are expected to be of the order of 10"" cm' s"', lead- 
ing to critical densities of 10' to 10* cm'. 
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3.1. Can we constrain the size of the CH2DOH 
emission ? 

We plotted the rotational diagrams of CH2DOH towards 
IRAS4A (the source with the largest number of detected lines) 
assuming different sizes for the source emission: 10", 15", 20", 
25" and 30". The data have been acquired at frequencies span- 
ning 89 GHz to 223 GHz, implying beam sizes between 30" 
and 11" (Table|6j. The different hypothesis on the source size 
will then correct differentially for the dilution depending on the 
frequency of the transition: if the beam size is smaller than the 
source, no correction is done, whereas dilution has to be cor- 
rected for observations with beam sizes bigger than the source. 
Methanol has the specificity to show no monotonic relation 
between the upper energy and the frequency of the transition 
(contrarily to CO for example), and thus the dilution correction 
has quite an unpredictable effect on the rotational diagram. 

For each assumed source size we calculated the reduced 
of the linear fit in the rotational diagram. If the assumed 
source size is the only origin of the scattering in the diagram, 
this reduced is minimum when the assumed source size is 
close to the real source size. 

The results of this study are presented in Table |2] The re- 
duced decreases when increasing the source size, suggesting 
that the source is extended. Nevertheless, this whole study lies 
on three assumptions: the emission is homogeneous on the ex- 
tent of the source, the source size is the same for all transitions, 
and scattering in the rotational diagram is only caused by the di- 
lution effect. These three hypothesis may not be valid. Indeed, 
the different transitions may be emitted by different regions, for 
instance the high energy transitions may originate in warmer 
and less extended regions than low lying transitions. The scat- 
tering in the diagram may also be caused by opacity or non- 
LTE effects. It may thus be unrealistic to derive the source size 
by this method. Only interferometric observations may help to 
solve this issue. 

Nevertheless, it is worth studying the uncertainties led by 
the source size assumption on the column density of the various 
isotopomers. Table|2]presents the column densities of methanol 
isotopomers versus the source size, obtained as follows. For 
CH3OH, column densities can only be estimated accurately for 
sizes of emission smaller than 15", because no observation was 
made with a beam larger than 15". CH3OD and CHD2OH col- 
umn densities were determined using the rotational tempera- 
ture derived for CH2DOH, for which we have the largest num- 
ber of transitions. Column densities change by a factor 2 to 3 
according to the size of the source. On the contrary, the ratios 
between isotopomers vary only slightly with the source size 
(cf. Table|5Jl. In the following, we will thus present the column 
densities assuming a source size of 10". Finally, care should 
be taken in the use of the fractionation ratios relative to the 
main isotopomers (CH3OH or H2CO), as these species might 
be optically thick. On the contrary, ratios between deuterated 
isotopomers are likely not to be affected by such problem and 
are thus reliable. 




10 20 30 40 

Upper energy (cm"^) 



Fig. 1. Rotational diagrams for deuterated formaldehyde for the 
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Table 2. Column densities for methanol isotopomers towards IRAS4A, derived by fixing the rotational temperature to the one 
derived for CH2DOH, and ratios between isotopomers. 



Source 


T,-o, 


CH2DOH 


2 

■^red 


CH3OD 


CHD2OH 


CH3OD 


CHD2OH 


size 


(K) 


(xlO'^ cm--) 




(xlO'^^cm-2) 


(xl0"cm-2) 


CH2DOH 


CH2DOH 



10" 


27.1±1.3 


4.3±04 


13.5 


3.1+0.7 


11+1.7 


(7.2+1. 8)xlO- 


-2 


0.26+0.05 


15" 


31.3±1.7 


2.4±0.3 


8.5 


1.6+04 


6.3+1.0 


(6.7+1.9)xl0- 


-2 


0.26+0.05 


20" 


36.2±2.3 


1.8±0.2 


5.2 


1.2+0.3 


5.0+0.9 


(6.7+1.8)xl0- 


-2 


0.28+0.06 


25" 


42.2+3.2 


1.7+0.2 


3.2 


1.4+0.3 


4.4+0.8 


(8.2+2.0)xl0- 


-2 


0.26+0.06 


30" 


46.2±3.8 


1.7+0.2 


2.9 


1.5+04 


3.9+0.7 


(8.8+2.6)xl0- 


-2 


0.23+0.05 



Table 3. H2CO, HDCO and D2CO rotational temperatures and column densities, assuming a 10" source. 



Source 




H2CO" 




HDCO 




D2CO 




T.-0, 




T.-0, 


N,„, 




T.-0, 


N,o, 




(K) 


(xlO'-'cm--) 


(K) 


(xlO'-'cm^-) 


(xlO'-'cm^^) 


(K) 


(xlO'^cm^^) 


IRAS4A 


20+2 


1.7+04 


7.3+0.3 


3.2+0.6 


7.2+0.2 


5.0+0.3 


19+4 


IRAS4B 


38+6 


1.7+0.5 


8.2+0.5 


1.9+0.4 


1.2±0.3 


7.9+0.5 


7.0+1.5 


IRAS2 


20+2 


0.84+0.22 


9.3+0.6 


1.3+0.3 


0.82±0.13 


19.6+3.6 


4.0+1.5 


L1448N 


14+1 


1.1+0.2 


7.9+0.5 


1.0+0.3 


0.46±0.057 


4.9+0.8 


8.2+4.4 


L 1448mm 


16+1 


0.42+0.11 


8.7+0.6 


1.1+0.3 


0.57±0.063 


5.7+0.3 


9.1+2.3 


LI 157mm 


12+1 


0.20+0.06 


10.0+0.9 


0.29+0.09 


0.24±0.034 




< 1.5 


L1527 


12+1 


0.42+0.17 


4.8+0.3 


6.0+1.8 


1.1±0.16 


5.1+0.3 


14.9+4.2 



"Reanalysis of the data from lMaret et alJ )2004) with the rotational diagram method {cf. text). ''Assuming the same rotational temperature as 
H2CO. 

Table 4. Methanol rotational temperatures and column densities, assuming a 10" source size. 



Source 


CH3OH 


CH2DOH 


CH3OD 


CHD2OH 




Trot 


N,„, 


Trot 


N,„, 


N,„, 


N,o, 




(K) 


(xlO"* cm-2) 


(K) 


(xlO"* cm-2) 


(xl0"cm-2) 


(xlO"* cm-2) 


IRAS 4A 


38.0+3.2 


6.9+14 


27.1+1.3 


4.3+0.4 


3.1+0.7" 


1.1+0.17" 


IRAS 4B 


84.9+16.8 


8.0+2.7 


15.6+0.5 


2.9+0.2 


1.1+0.20'' 


0.90+0.10" 


IRAS 2 


207+48.1 


10.1+3.7 


55.3+4.9 


5.2+0.8 


< 8.0 


2.1+0.4" 


L1448N 


27.5+4.5 


1.2+0.47 


27.5+4.5" 


2.1+0.6* 


<8.0 


< 10 


L1448mm 


76.3+19.1 


1.6+0.70 


76.3+19.1" 


11.3+4.7* 


40+16* 


<6.3 


LI 157mm 


112+201 


1.9+5.2 


112+201" 


10.3+27.7* 


22+66* 


< 9.6 



" Temperature fixed to the rotational temperature of CH3OH. 

* Should be taken as an upper limit, as the rotational temperature is likely to be overestimated. 
" Temperature fixed to the rotational temperature of CH2DOH. 



3.2. Formaldehyde and methanol fractionation 

Figure[l]presents the rotational diagrams derived for the deuter- 
ated isotopomers of formaldehyde towards the seven sources 
of the sample assuming 10" source sizes. Table |3 lists the de- 
rived rotational temperatures and column densities. For con- 
sistency with the analysis of HDCO and D2CO, the H2CO 
column density has been recom puted with t he rot ational di- 
agram method from the data of iMaret et alJ (|2004), with the 
ortho/para ra tio fixed to its statistical value (contrarily to the 
study of Maret et alJ ( Eo04.) where this ratio was considered as 



a free parameter), and taking into account the line opacities 
for the transitions where H2''^CO has been observed. Note that 
the column densities inferred by treating the ortho/para ratio 
as a free parameter are less than a factor 2 different from our 
values ( Maret et alj2 004). The derived rotational temperatures 
for HDCO and D2CO are in some cases different from the one 
derived for H2CO, but it should be noticed that the transitions 
only span a small interval of upper energies and thus do not 
constrain very well the rotational temperature. In order to get 
a sense of the uncertainty implied by these different rotational 
temperatures, we also computed the HDCO column densities 
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assuming the same rotational temperature as H2CO (cf Table 
Results were found to be in most cases within a factor of 
two uncertainties, and the fractionation ratios mostly appear to 
be consistent within the error bars (cf Table|5] footnote a). 

Fig. |2] presents the rotational diagrams of the deuterated 
forms of methanol for the three sources IRAS4a, IRAS4b 
and IRAS2 respectively, again assuming a source size of 10". 
The CH3OD and CHD2OH column densities were computed 
assuming the same rotational temperature as CH2DOH, for 
which many more lines were observed in all sources. As for 
form aldehyde, the CH^ OH column density was recomputed 
from lMaret et alJ ll20Q5h using rotational diagrams. 



4. Discussion 

In this section, we discuss the results of the observation of 
deuterated formaldehyde and methanol towards the sample 
of low-mass class protostars. Table |5l summarizes the ob- 
served deuterium fractionations. For com parison, we added the 
fractionations me a sured i n IRAS 16293 dLoinard et al.l 
IParise et all [20021 12004 . iRoberts et alJ (|2002l) observed the 
HDCO fractionation in a number of protostellar cores, in- 
cluding L1448mm and L1527, using the Kitt Peak 12 m tele- 
scope. They found values sensibly lower than ours in these 
two sources (resp. 0.069 and 0.066). This may come from the 
fact that the 12 m telescope beam size is 2.5 times bigger than 
the 30 m, so that their observations encompass the surrounding 
cloud where the fractionation is expected to be lower 

The main result of this observational study is the discovery 
of a high fractionation ratio for formaldehyde and methanol for 
each source where these molecules were detected. The mea- 
sured fractionation ratios are indeed similar to those observed 
in IRAS 16293. In the sources where the isotopes were not de- 
tected, the derived upper limits do not exclude a similar deuter- 
ation. The present observations confirm that IRAS 16293 is not 
an exception regarding deuteration. Other studies also confirm 
that its chemistry is not peculiar Indeed Bottinelli et al. ( 2004) 
observed complex molecules in the Hot Core of IRAS4A, 
molecules that were detected for the first time in a low-mass 
protostar towards IRAS 16293 jCazaux et alJ2003h . 

We explore in the following paragraph the correlations be- 
tween the deuteration and other parameters. We then compare 
the observed fractionations to the fractionations predicted by 
grain chemistry models. 

4.1. Correlations 

Figure |3l shows the measured D2CO/H2CO and 
CH2DOH/CH3OH as functi on of the bolom etric luminos- 
ity and the Ls,n,„/L[,oi (from IXiidre et al.l '2000). The figures 
show no obvious correlations between the measured frac- 
tionation and these two quantities, despite a variation of 
one order of magnitude in the D2CO/H2CO ratio. This is 
indeed consistent with the hypothesis that the deuteration is 
a memory of precollapse phase, rather than a present day 
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Fig. 2. Rotational diagrams for deuterated methanol in 
IRAS4A, IRAS4B and IRAS2 (for the sake of clarity the curves 
for CH3OD and CHD2OH have been tianslated by -2 and -7 
along the y-axis). Error bars correspond to the error bars on 
the flux as stated in Tab. |6l and |2l Solid lines (CH2DOH) cor- 
respond to fits to the data, with two free parameters (Trot and 
Ntot). Dashed lines correspond to fits with only Ntot as a free 
parameter, and Tiot set to the value derived for CH2DOH (cf 
text). 



product ("e.g. ICeccarelli etallbOOlh . If one cannot therefore 
use the fractionation to distinguish the evolutionary status of 
the sources, it is likely that the fractionation gives insights of 
the precollapse phase of each source. 

Figure |4] presents the D2CO/H2CO fractionation ratio (as 
well as the D2CO/HDCO ratio in order to get rid of any poten- 
tial opacity eff'ect) versus the CO depletion in the sources of our 
sample. The CO depletion was computed from the CO densities 
in the outer envelope derived bv lj0rgensen et al. ( 2002 ), an d as- 
suming a canonical CO abundance of 8.4x10"'^ iFrerking et alJ 
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Table 5. Deuterium fractionation for formaldehyde and methanol, for an assumed source size of 10". Error bars are associated 
with a 90% confidence interval (see appendix). 



Source 


HDCO" 


D2CO 


D2CO 


CH2DOH 


CH3OD 


CHD2OH 


CHD20H 




HjCO 


H2CO 


HDCO 


CH3OH 


CH3OH 


CH3OH 


CH2DOH 


IRAS4A 


20+" ' ^^ 

"•^"-0.09 


12+°°* 


62+° " 

-0.26 


65+°-^^° 
"•"•^-0.21 


047+° °^'' 
"•"^'-0.021 


17+°°'' 

^ ' -0.06 


26+° °* 

"•^"-0.07 


IRAS4B 


13+°'° 

" -^-0.06 


046+° °^'* 


39+°-^^ 


43+°-38 
^ -^-0.20 


016+° °"' 


13+°'- 


31+°°' 

"•■^^-0.06 


IRAS2 


17+°'^ 


052+° °'"* 


33+°-2* 


62+°-^' 
-0.33 


< 0.08 


25+°-2'' 


41+°" 


L1448N 


094+0.062 

"■"^^-0.050 


077+° °" 

' ' -0.069 


93+'-"' 


< 1.8 


< 0.67 


< 8.3 




LI 448mm 


29+°-2^^ 


24+° '** 

"•^^-0.12 


QQ^+0.73 


< 7.1 


<2.5 


< 3.9 




LI 157mm 


16+°" 


< 0.08 


< 0.52 


< 5.4 


< 1.2 


< 5.1 




L1527 


1 7+2'' 


4.4.+0.60 
"•^^-0.29 


0.28!°:^^ 










IRAS 16293 


0.15±0.07 


0.05±0.025 


0.3±0.2 


bn 17+0.38 
"•■-"-0.19 


*o.oi8!°:°?^ 


''O 074+° °'*"' 
"•"'^-0.044 


*Q 21+0.11 

"•^^-0.10 



°HDCO/H2CO ratios obtained when assuming the same rotational temperature for HDCO as for H2CO are resp. 0.075+° °;j^, 0.080+° °5*, 
0.11!°°^, 0.043!°°2i, 0.15!°i°, 0.14+°^', 0.33!^^^^ 

HDCO/H2CO ratios obtained assuming a 15" source are resp. 0.14;^°;^^, 0.07+°;°^, 0.14;^°;^^, 0.08!°;°^, 0.23!°;[^, OA^lf^, 0.89!,!,;^^. They are 

thus consistent with the ratio s derived for a 10" source, within the error bars. 

* Error bars recomputed from lParise et alJ^2004^ for a 90% confidence interval (see appendix). 
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Fig. 3. Fractionation of formaldehyde and methanol as a function of bolometric luminosity and the Ls,„mfLboi (from lAndre et alJ 

be 



Il982h . We also p lotted in this figure th e result for prestellar 
cores observed bv lBacmann et alJ ( l2002l) . 

The first thing to notice is that Class sources and prestel- 
lar cores with similar CO depletion factors do present also sim- 
ilar D2CO/H2CO ratios. This is a further suggestion that the 
deuteration of formaldehyde in Class sources occurred dur- 
ing the precollapse phase, and/or that the physical conditions 
in the outer envelopes of Class sources are very similar to 



prestellar cores. On the other hand, and surprising enough, the 
two sources with the largest D2CO/H2CO ratio are those with 
the lowest CO dep letion factor. Indeed chemical models (e.g. 
iRoberts & Millaj EoOO) predict that formaldehyde fractiona- 
tion should increase with CO depletion. However, we inquire 
in the following if this could be caused by the optical thick- 
ness of H2CO, leading to artificially high D2CO/H2CO ratios. 
Both sources (L1527 and L1448mm) present evidences of very 
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Fig. 4. Doubly-deuteration of formaldehyde versus CO deple- 
tion. Class sources (this work) appear as diamonds and 
prestellar cores ( Bacmann et al. 2002|) as triangles. The re- 
versed triangle is the upper limit for LI 157mm (this work). 



optic ally thick CO emissio n. A study of C'^O and C'^O emis- 
sion ( |j0rgensen et al.l2002h towards a sample including all our 
sources showed that L1527 and L 1448m m are the sources fo r 
which C'^O is the most optically thick (ileirgensen et alJ2 0Q2). 
H2CO is thus more likely to be optically thick (as suggested 
by the correlation between the H?CO and CO abundances, 
iMaret et all l2004 . iMaret et all ( |2004 have observed some 
transitions from Hj^'CO in a few sources. Unfortunately, no line 
was observed for LI 527, and only an upper limit was derived 
for one line on L1448mm (r < 2). In the case of L1448mm, we 
thus might have underestimated the H2CO column density by 
a factor as high as ~ 2.3. For L1527, the opacity is likely 
to be even large r as su ggested by the C'^O/C'^O measured by 
|j0rgensen et alJ ( l2002h . However, we conclude from this study 
that this opacity effect cannot fully account for the discrepancy 
of 2 orders of magnitude observed in D2CO/H2CO compared to 
prestellar cores. Indeed, the D2CO/HDCO, which is supposed 
to be less affected by opacity effects than the ratio with the main 
isotopomer, is still too high for one of the sources (L1448mm). 
Hence, the conclusion is while there seems to be a correlation 
between CO depletion and fractionation in prestellar cores, no 
such correlation is apparent in our sample of YSO (no corre- 
lation is found either between the fractionation and the N2H^ 
abundance, tracing the N2 depletion). This difference may im- 
ply that the CO depletion in the circumstellar envelopes is af- 
fected by outgassing due to the heating by the newly formed 
star. 

Finally, our analysis, based on rotational diagrams, allows 
only to estimate the fractionation averaged on the beam. It will 
be important to understand what the contribution of the dif- 



ferent regions (warm envelope, cold envelope, etc) is. New in- 
sights in this respect may only be possible by interferometric 
observations, to observe the relative fractionation of formalde- 
hyde and methanol in the two regions ; hot core and outer en- 
velope. 

4.2. Comparison to grain chemistry models 

Grain surface chemistry is thought to be responsible for the 
formation of the abundant methanol observed in hot cores, for 
two main reasons. First, gas-phase chemistry models fail to re- 
produce the CH^OH abundance by 2 to 3 orders of magnitude 
in Orion (Menten et al. 1988). Recent investigations have con- 
firmed that no gas-phase process can efficiently form methanol 
llHerbstlbQOsI) . Second, methanol is one of the most abundant 
species observed in the grain mantles (after water) , with abun- 
dance s as high as 30% of the water abundance ( IDartois et alJ 
119901) . 

Methanol is thought to be formed on the grains by suc- 
cessive hydrogenations of CO. An intermediate product of 
these reactions is formaldehyde. Deuterium fractionation stud- 
ies of formaldehyde and methanol thus provide a useful tool 
for confirming the grain chemistry scenario. To understand if 
formaldehyde and methanol are formed simultanuously on the 
grains, and to unveil a possible contribution of gas-phase pro- 
cesses in the formation of formaldehyde, we compared the ob- 
served fractionations to the predictions of the grain chemistry 
model from Stantcheva & Herbst (2003). This model is based 
on the resolution of the rate equations governing the abundance 
of each species on the grain, and uses the master equation ap- 
proach in the cases where only a small number of molecules 
are present on the surface. 

Fig. |5l presents, in dashed lines, the HDCO, D2CO, 
CH2DOH, CH3OD et CHD2OH fractionations predicted by the 
grain model, as a function of the gas-phase atomic D/H ratio at 
the time of mantle formation. The observed fractionations with 
their error bars have been superimposed for each source. This 
allows to infer the required D/H ratio required for the formation 
of each molecule. 

For a compa rison. Fig. |5] also pre s ents the same study 
for IRAS16293 jCeccarelli et alJll998t lT,oinard et alj 12000: 
IParise et alJ2004 '). 

A comparison between formaldehyde on the one hand and 
methanol on the other hand must be led with care as it is very 
dependent on the possible opacity of the main isotopomers 
H2CO and CH3OH. On the opposite, a comparison between 
different deuterated isotopomers of methanol should not suffer 
from such problem. We attribute to this opacity issue the fact 
that, depending on the source, either HDCO or D2CO is com- 
patible with methanol fractionation. 

Several conclusions can be drawn from this comparison : 

• Whenever they are correctly constrained, the CH2DOH 
and CHD2OH fractionations are consistent, i.e. they have a 
common range of possible atomic D/H ratio (IRAS 16293, 
IRAS4A, IRAS4B et I RAS2). Tha t had already been ob- 
served for IRAS 16293 dParise etalJEoOZ 20oA for which 




Fig. 5. Comparison between the observed fractionations (triangles are for formaldehedyde species whereas diamonds are for 
methanol species) towards the seven sources of our sample (sohd lines) and the predictions of the grain model (dashed hnes, 
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CD3OH fractionation was also consistent with the latter two 
isotopomers. The present new observations thus add further 
support to the hypothesis that methanol is formed on the grain 
surfaces. 

• Whenever it is contrained, CH3OD appears to be un- 
derabundant with respect to the other methanol isotopomers. 
Th is had alrea d y bee n noticed for the case of IRAS 16293 
by 'Paris e et all (|2.002"), who concluded that CH3OD may be 
destroyed in the gas phase by protonation more rapidly than 
the other isotopomers, owing to the strong electronegativity 
of the oxygen atom. This conclusion was also suggested by 
the theoretical study of Osarnura et al. (2 004) . However, some 
very re cent laboratory experim ents suggest another possibility. 
Indeed, 'Naga oka et alJ ( l2005l) show that deuterated methanol 
species with an OD bond are very inefficiently formed on ice 
surfaces when exposing CH3OH to H and D atoms, while 
isotopomers with CD bonds are abundantly formed. This might 
also explai n the non-detection of HDO in ices of p rotostellar 
envelopes dOartoisetalJ l2003t IPariseetalJ l2003l) and the 
fractionation of water observed to be one order of magnitude 
lower than CH2DOH/CH3OH in the hot corino of IRAS 16293 
in agr eement with these solid state observations dParise et alJ 

• HDCO appears to be less abundant than D2CO with 
respect to the model predictions (except for L 11 57mm, but 
for this source, the error bars are important). No prefer- 
ential destruction in the gas phase could explain this dis- 
crepancy (as it was the case for CH3OD). This might be 
caused by H-D substitution via hydrogen abstraction reactions 
(HDCO -H H ^ DCO + H2), which are not included in the con- 
sidered chemistry model. These reactions would favour the for- 
mation of D2CO by abstra cting an hydrogen to HDCO and 
deuterating DCO llTielensll9 83). 

Alternatively, this might be a hint that formaldehyde forma- 
tion and fractionation is not dominated by grain chemistry pro- 
cesses. We will discuss more thoroughly the gas-phase forma- 
tion and fractionation of formaldehyde in the next section. 

4.3. A gas-phase formation of formaldehyde ? 

As discussed in the previous section, there is now good ev- 
idence that methanol is formed on the grains, and the data 
of this survey add further support to this hypothesis (as pre- 
sented in the previous paragraph). Yet, the situation is not so 
clear for formaldehyde. After discussing the spatial origin of 
the HDCO and D2CO emission, we examine the implications 
of the present observations on the formaldehyde formation pro- 
cess. 

Contrarily to methanol that shows quite high rotational 
temperatures, and thus may be arising in the hot corinos, 
deuterated formaldehyde shows low rotational temperatures, 
rather pointing to a n origin in the cold en velope of the pro- 
tostars. The study bv lCeccarelli et all 1I2OOII) of the distribution 
of the D2CO 4o,4-3o,3 line (E„p=19.4cm~\ a rather low energy 
line) towards IRAS 16293-2422 shows that it is extended over 



the entire envelope, and that the D2CO/H2CO ratio is also large 
across the envelope. These authors also noticed that while the 
D2CO emission is well correlated with the continuum emis- 
sion tracing the envelope, it does not show any enhancement 
in correspondence of the shocks sites in the region caused by 
the interface between the outflow and the envelope itself. They 
concluded that D2CO is abundant in the cold outer envelope, 
where the dust temperature is lower than 100 K. Of course, this 
does not exclude the possibility of abundant D2CO also in the 
hot corino region, but the available observations were unable 
to detect it. Note that Fuente et al. (2005) found, from infer- 
ferometric observations, abundant D2CO in the hot core of the 
intermediate mass protostar NGC7129-F1RS2. Similarly to the 
case of IRAS 16293-2422, the present observations are unable 
to disentangle a possible contribution from the hot corino, as 
well as from the shocks around these sources. But they cannot 
exclude these contributions either Given the uncertainty on the 
location of the D2CO, in the following we analyze both possi- 
ble situations, hot corino versus cold envelope. 

T he high observed deuterium fractionation of formalde- 
hyde (jCecc^emetanillll |Loi^^ and this pa- 
per), if arising in the hot corino, excludes an efficient produc- 
tion of formaldehyde in the warm gas, for which a much lower 
deuteration would be expected. Hence these high formaldehyde 
abundance and fractionation must reflect the cold phase pre- 
ceeding the hot core phase accompanying the formation of the 
protostar Gas phase models typically produce formaldehyde 
abundance which range from about 10"'' at early times to 10"** 
at late times (e.g. in steady state). Observed abundance in cold 
dark clouds are typically 10"** and in good agreement with the 
latter value. However, abundances are observed t o be higher 
in the hot co rino around protostars (~ 10"^- 10'^. iMaret et alJ 
l2004 12005I) . Hence, models based upon gas phase produc- 
tion of the formaldehyde in hot cores have to selectively de- 
plete the formaldehyde in order to keep the abundance in ffie 
dark phase low. While all species are expected to freeze out 
equally efficiently at the low temperature of dust grains in dark 
clouds, perhaps, CO is returned to the gas phase by evapora- 
tion but formaldehyde is left behind. Then, in order to increase 
the formaldehyde abundance by a factor of 100, the CO has to 
cycle a hundred times between the gas phase and the ices. This 
seems a bit contrived. The only gas phase possibility seems 
therefore formation at early times, when H2CO abundance is 
10"^, and immediate depletion and cold storage in unreactive 
ices. This seems to be in contradiction with the depletion stud- 
ies of dark cloud cores (Tafalla et al. 2002). Furthermore, the 
bulk of the formaldehyde would be formed under undepleted 
conditions and hence the D2CO and HDCO fractionation is dif- 
ficult to understand (Roberts et al. 2003). 

We now study the implications of our observations if the 
HDCO and D2CO mainly originate in the cold envelope. In 
this case, we must inquire whether gas phase processes at 
present times can account for their formation. The discrepancy 
of our fractionation observations for formaldehyde with the 
grain chemistry models (namely the fact that HDCO seems to 
be steadily underabundant compared to D2CO) might indeed 
reveal that the formation of formaldehyde is not dominated 
by H and D additions to CO on grain surfaces. In this cold 
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gas, the absolute H2CO abundance is well explained by 
gas-phase models. Models of formaldehyde fractionation in 
prestellar cores -whose chemical and physical conditions are 
similar to those in the outer cold envelopes of Class sources 
jMaret etalJi2004l: |j0rgensen et al I I2OO5I)- seem approxima - 
tively to reproduce the observed ratios ("Robe rts et al]l2004) . 
However, the predictions depend on some key reactions, whose 
rates, products, and branching ratios are relatively uncertain, 
and recent studies seem to show that gas-phase models may 
have more difficulties to reproduce the high deuteration than 
expected. Indeed, lOsamura et alJ fcoOSh recently presented 
quantum chemical calculations of the protonation/deuteration 
of H2CO and deuterated isotopologues. These computations 
show that the deuteration of H2CO by and deuterated 
isotopomers will mainly happen on the oxygen end of the 
molecule, due to the higher electronegativity of oxygen. It is 
generally expected that dissociative electron recombination of 
H2COD+ will not lead to the formation of HDCO, so that the 
fractionation due to reactions of H2CO with H2D^ and other 
isotopomers is quite inefficient. However, given the recent 
surprising experimental results on the product distribution in 
the dissociative reco mbination reaction of N2H^ and CH3OH2 
jGeppert & alJ l2006l) . an experimental verification of the 
product distribution of this key reaction (H2COD^ -1- e ) would 
be very valuable. Another remaining possibility in this context 
is that H2CO is formed in the gas-phase and then frozen out 
on dust grains where it can be later fractionated by abstraction 
reactions on the grain surfaces. Some efficient desorption 
processes are then required to release the molecules back to 
the gas-phase. 

In summary, pure gas-phase models seem to face problems 
to explain both formaldehyde absolute abundance and fraction- 
ation at the same time. Unless the dissociative recombination of 
H2COD+ leads to the formation of HDCO, H-D substitutions 
on grains seem to be required to account for the high fraction- 
ation if formaldehyde is formed by gas phase processes. 

5. Conclusion 

We have presented observations of deuterated formaldehyde 
and methanol towards a sample of low-mass protostars. We de- 
tected HDCO, D2CO and CH2DOH in all sources. However, 
CH3OD and CHD2OH were only detected towards 2 and 3 
sources respectively. We analyzed the data using population 
diagrams. These observations show that IRAS 16293-2422 is 
not an exception concerning deuteration. The fractionations de- 
rived towards the source of our sample are indeed similar to 
those observed in IRAS 16293. 

These observations are useful to pin down the contribu- 
tions of grain surface and gas phase chemistry for the forma- 
tion of formadehyde and methanol. It appears that methanol 
fractionation is consistent with grain surface schemes, whereas 
formaldehyde formation cannot be fully constrained from these 
data. Two possibilities arise, either that formaldehyde forma- 
tion is dominated by gas-phase reactions, later fractionated on 
the grains through abstraction reactions, or that formaldehyde 



is formed on the grains, with abstraction reactions altering the 
expected D2CO/HDCO ratio. 

This study is limited by different factors. All the transitions 
were observed with different sized beams, with a best spatial 
resolution of 10". It is in this case very difficult to derive the 
abundance of the molecules if one does not know the source 
size. Moreover, the spatial resolution is too poor to derive the 
fractionation in the different regions of the sources (warm in- 
ner envelope, cold outer envelope, outflows ...). Interferometry 
is the only chance to go further in this study. In this respect, 
ALMA will enable considerable progress. 

Acknowledgements. B. Parise wishes to thank the IRAM 30 m opera- 
tor Manuel Ruiz for his precious help during the observations. 
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CH2DOH/CH3OH 

Fig. 6. Repartition function for the CH2DOH/CH3OH ratio 
(solid line). Overlaid in dash-dot hue is the repartition func- 
tion computed assuming that the ratio is Gaussian, with a width 
computed by the standard error propagation. Vertical lines de- 
limitate the 90% confidence intervals for the two distributions. 



Appendix : Computation of error bars for the frac- 
tionation ratios 

The ratio of two Gaussian distributed variables X and Y is 
Gaussian only when crx « X and cry « Y. When the errors 
are too large, deviations from Gaussianity are observed and it 
is thus necessary to compute accurately the repartition func- 
tion of the random variable R=X/Y. This repartition function is 
given by: 

fx(uy)fY(y) lyldy 

DO 

This distribution is generally asymmetric, and thus leads to 
asymmetric confidence intervals. We chose in this paper to give 
intervals associated to a 90% confidence level. As an example, 
we present here the comparison with the repartition function for 
the CH2DOH/CH3OH ratio in IRAS2, assuming that CH2DOH 
and CH3OH have a Gaussian distribution. 
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ONLINE MATERIAL 



L1527 



Frequency 


Fup 




Tmb 


Av 


GHz 


cm"' 


K.km.s-1 


K 


km.s"' 


HDCO 


134.285 


12.3 


0.75+0.08 


1.1 


0.7+0.1 


201.341 


19.0 


0.70±0.08 


0.98 


0.7±0.1 


256.585 


21.4 


0.77±0.10 


0.18 


0.9±0.1 


D2CO 


110.838 


9.3 


0.133±0.016 


0.243 


0.5+0.1 


166.103 


14.8 


0.219±0.038 


0.327 


0.6+0.1 


231.410 


19.4 


0.367±0.047 


0.601 


0.6+0.1 


276.060 


31.4 


<0.23 







Table 8. Main-beam intensities, peak temperatures and width 
for the observed transitions towards L1527. 



Table 9. Main-beam temperature upper limits at 3cr level for 
the deuterated methanol transitions observed towards L1448N, 
L1448mm and LI 157nam. 





L1448N 


LI 448mm 


LI 157mm 


Freq. 


/ Tmbdv 


/ Tmbdv 


/ T„,bdv 


GHz 


K.km.s-' 


K.km.s~' 


K.km.s-' 


CH2DOH 


89.251 


< 0.015 


<0.03 


<0.02 


89.275 


< 0.015 


<0.03 


<0.02 


110.105 


<0.03 


<0.05 


<0.04 


207.781 


<0.08 


<0.08 


<0.05 


223.071 


< 0.05 


<0.08 


<0.05 


223.107 


<0.05 


<0.08 


<0.05 


223.128 


<0.05 


<0.08 


<0.05 


223.131 


<0.05 


<0.08 


< 0.05 


223.131 


<0.05 


<0.08 


< 0.05 


223.154 


<0.05 


<0.08 


< 0.05 


223.154 


<0.05 


<0.08 


<0.05 


223.315 


<0.04 


<0.06 


<0.05 


223.422 


<0.04 


<0.06 


<0.04 


CH3OD 


110.263 


<0.03 


<0.05 


<0.03 


110.476 


< 0.02 


< 0.04 


< 0.03 


223.304 


< 0.04 


< 0.06 


< 0.05 


CHD2OH 


207.771 


<0.08 


<0.08 


< 0.04 


207.827 


<0.08 


<0.08 


<0.04 


207.864 


<0.08 


<0.08 


<0.04 


207.868 


<0.08 


<0.08 


<0.04 


207.869 


<0.08 


<0.08 


<0.04 
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NGC1333-IRAS4A NGC1333-IRAS4B NGC1333-IRAS2 




5 10 

Velocity (km/s) 



5 10 
Velocity (km/s) 



5 10 
Velocity (km/s) 



Fig. 7. HDCO lines for IRAS4a, IRAS4b, IRAS2. 
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Fig. 8. HDCO lines for L1448N, L1448mm, LI 157mm 
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NGC1333-IRAS4A NGC1333-IRAS4B NGC1333-IRAS2 
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Fig. 9. D2CO lines for IRAS4a, IRAS4b, IRAS2. 
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Fig. 10. D2CO lines for L1448N, L1448mm, LI 157mm 
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NGC1333-IRAS4A NGC1333-IRAS4B NGC1333-IRAS2 
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Fig. 11. CH2DOH lines for IRAS4a, IRAS4b and IRAS2. 
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Fig. 12. CH2DOH lines for IRAS4a, IRAS4b, IRAS2. 
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Table 6. Main-beam intensities, peak temperatures and width for the observed transitions towards NGC1333-IRAS4A, IRAS4B 
and IRAS2. The errors on the fluxes were computed as the quadratic sum of the statistical error and the calibration uncertainty 
(see text). Upper hmits are 3cr (see text). 









IRAS4A 






IRAS4B 






IRAS2 




Freq. 




/ Trabdv 


Tmb 


Av 


/ Tmbdv 


Tmb 


Av 


/Tmbdv 


Tmb 


Av 


GHz 


cm ' 


K.km.s ' 


K 


km.s ' 


K.km.s ' 


K 


km.s"' 


K.km.s-' 


K 


km.s"' 



HDCO 



134.285 


12.3 


0.95+0.10 


0.58 


1.5+0.1 


0.59+0.06 


0.43 


1.3+0.1 


0.43+0.04 


0.33 


1.2+0.1 


201.341 


19.0 


1.09±0.11 


0.73 


1.5+0.1 


0.70+0.07 


0.47 


1.4+0.1 


0.69+0.07 


0.53 


1.2+0.1 


246.924 


26.1 


0.98±0.11 


0.65 


1.6+0.1 


0.74+0.10 


0.57 


1.2+0.1 


0.50+0.08 


0.33 


1.4+0.2 


256.585 


21.4 


1.19±0.13 


0.94 


1.2+0.1 


1.38+0.15 


0.83 


1.6+0.1 


1.11+0.14 


0.92 


1.1+0.1 


268.292 


27.9 


0.82±0.13 


0.64 


1.2+0.1 


0.74+0.12 


0.44 


1.6+0.1 


1.00+0.16 


0.76 


1.1+0.1 


D2CO 


110.838 


9.3 


0.16±0.01 


0.12 


1.2+0.1 


0.095+0.011 


0.12 


0.8+0.1 


0.077+0.017 


0.062 


1.2+0.4 


166.103 


14.8 


0.33±0.04 


0.24 


1.3+0.2 


0.24+0.04 


0.22 


1.0+0.2 


0.095+0.022 


0.21 


0.4+0.1 


231410 


194 


0.37±0.05 


0.28 


1.3+0.1 


0.47+0.06 


0.41 


1.1+0.1 


0.63+0.07 


0.76 


0.8+0.1 


275.671 


20.9 








<0.16 






<0.26 






276.060 


31.4 








0.10+0.03 


0.46 


0.2+0.1 


0.44+0.10 


0.45 


0.9+0.2 


CH2DOH 


89.251 


17.2 


0.042±0.008 


0.014 


2.8+0.6 


0.011+0.005 


0.008 


1.3+0.8 


0.021+0.008 


0.006 


3.4+1.2 


89.275 


13.8 


0.033±0.007 


0.015 


2.0+0.5 


0.023+0.004 


0.015 


1.4+0.3 


< 0.013 






89.408 


4.5 


0.098±0.008 


0.054 


1.7+0.1 


0.113+0.008 


0.072 


1.5+0.1 


0.038+0.005 


0.015 


2.3+0.4 


110.105 


834 


0.054±0.014 


0.022 


2.3+0.6 


< 0.027 






0.050+0.009 


0.025 


1.8+34 


133.847 


18.3 


0.084±0.022 


0.046 


1.7+0.4 


< 0.077 






0.090+0.017 


0.045 


1.9+0.3 


133.873 


21.7 


0.064±0.019 


0.025 


2.4+0.7 


0.035+0.009 


0.030 


1.1+0.2 


0.104+0.020 


0.034 


2.9+0.5 


133.882 


33.6 


0.054±0.019 


0.018 


2.6+0.9 


< 0.039 






0.048+0.026 


0.016 


2.9+2.0 


133.893 


27.4 


< 0.052 






0.055+0.018 


0.024 


2.1+0.7 


<0.10 






133.930 


27.4 


< 0.052 






<0.039 






0.080+0.022 


0.025 


3.0+0.7 


134.066 


8.9 


0.157±0.022 


0.107 


1.4+0.2 


0.197+0.024 


0.160 


1.2+0.1 


0.115+0.021 


0.042 


2.6+0.4 


134.112 


20.2 


0.068±0.015 


0.071 


0.9+0.2 


< 0.052 






< 0.065 






207.781 


15.9 


0.094±0.017 


0.067 


1.3+0.2 


0.089+0.017 


0.024 


3.6+0.6 


0.144+0.029 


0.067 


2.0+0.5 


223.071 


33.6 


0.116±0.022 


0.059 


1.9+0.3 


0.042+0.011 


0.034 


1.1+0.3 


0.229+0.035 


0.081 


2.7+04 


223.107 


35.1 


< 0.10 






0.030+0.010 


0.024 


1.2+0.4 


0.234+0.041 


0.098 


2.3+04 


223.128 


40.8 


< 0.083 






0.046+0.013 


0.037 


1.2+0.4 


0.189+0.055 


0.073 


2.4+0.8 


223.131 


79.4 


0.015±0.01 


0.025 


0.6+0.2 


0.02+0.01 


0.02 


0.9+0.3 


0.04+0.03 


0.016 


2.3+1.2 


223.131 


79.4 


0.015±0.01 


0.025 


0.6+0.2 


0.02+0.01 


0.02 


0.9+0.3 


0.04+0.03 


0.016 


2.3+1.2 


223.154 


60.8 


0.046±0.02 


0.02 


2.0+0.5 


0.02+0.01 


0.015 


1.2+0.5 


0.13+0.03 


0.05 


2.6+0.4 


223.154 


60.8 


0.046±0.02 


0.02 


2.0+0.5 


0.02+0.01 


0.015 


1.2+0.5 


0.13+0.03 


0.05 


2.6+0.4 


223.315 


40.8 


0.061+0.016 


0.040 


1.4+0.4 


0.083+0.015 


0.032 


2.5+0.4 


0.192+0.032 


0.061 


3.0+0.4 


223.422 


33.6 


0.126±0.024 


0.061 


2.0+0.4 


0.062+0.010 


0.042 


1.4+0.2 


0.229+0.038 


0.064 


3.4+0.5 


CH3OD 


110.189 


7.8 


<0.03 






0.022+0.009 


0.054 


0.5+0.2 


< 0.025 






110.263 


10.8 


0.041+0.014 


0.018 


2.0+0.7 


<0.03 






<0.02 






110.476 


154 


<0.02 






<0.02 






<0.02 






133.925 


6.0 


0.076+0.020 


0.024 


2.7+0.8 


0.046+0.015 


0.038 


1.1+0.5 


<0.05 






223.309 


26.8 


<0.07 






0.035+0.009 


0.029 


1.2+0.3 


<0.15 






226.539 


22.7 


<0.06 






<0.04 






<0.17 






CHD2OH 


83.1292 


17.0 


<0.03 






0.018+0.005 


0.026 


0.7+0.2 


< 0.021 






83.2895 


4.2 


0.009+0.002 


0.016 


0.5+0.1 


0.024+0.004 


0.027 


0.9+0.2 


< 0.021 






83.3036 


10.3 


0.015+0.006 






0.011+0.004 


0.012 


0.9+0.3 


< 0.021 






207.771 


33.6 


0.066+0.040 


0.022 


2.8+1.7 


0.089+0.025 


0.025 


3.4+0.9 


0.115+0.041 


0.029 


3.7+14 


207.827 


42.6 


<0.07 






0.097+0.027 


0.021 


4.4+1.0 


0.087+0.045 


0.016 


5.2+2.3 


207.864 


68.4 


0.062+0.030 


0.036 


1.6+0.7 


< 0.035 






0.088+0.033 


0.035 


2.3+1.0 


207.868 


53.5 


0.12+0.04 


0.03 


3.7+0.8 


0.034+0.011 


0.02 


1.8+0.4 


0.13+0.03 


0.03 


4.8+0.9 


207.869 


53.5 


0.12+0.04 


0.03 


3.7+0.8 


0.034+0.01 1 


0.02 


1.8+0.4 


0.13+0.03 


0.03 


4.8+0.9 
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Table 7. Main-beam intensities, peak temperatures and width for the detected transitions towards L1448N, L1448mm and 
L 11 57mm. The errors on the fluxes were computed as the quadratic sum of the statistical error and the calibration uncertainty 
(see text). Upper limits are 3cr. Upper limits for lines not detected in any of these three sources can be found in Tab.|5] 



L1448N L1448mm L1157mm 



Freq. 




J Tmbdv 


Tmb 


Av 


J Tmbdv 


Tmb 


Av 


J Tmbdv 


Tmb 


Av 


GHz 


cm"' 


K.km.s-' 


K 


km.s"' 


K.km.s-' 


K 


km.s"' 


K.km.s-' 


K 


km.s"' 


HDCO 


134.285 


12.3 


0.32+0.04 


0.26 


1.1±0.1 


0.38+0.06 


0.38 


0.9+0.1 


0.11+0.01 


0.13 


0.8+0.1 


201.341 


19.0 


0.36±0.05 


0.38 


1.0±0.1 


0.53+0.06 


0.97 


0.6+0.1 


0.10+0.03 


0.12 


0.8+0.1 


246.924 


26.1 


0.35±0.04 


0.30 


1.1±0.1 


0.61+0.08 


0.63 


0.9+0.1 


0.26+0.03 


0.20 


1.2+0.1 


256.585 


214 


0.60±0.10 


0.47 


1.2±0.2 


0.79+0.10 


1.11 


0.7+0.1 


0.18+0.03 


0.18 


0.9+0.1 


268.292 


27.9 


0.34±0.08 


0.24 


1.4+0.3 


0.40+0.08 


0.66 


0.6+0.1 


0.15+0.05 


0.20 


0.7+0.4 


D2CO 


110.838 


9.3 


0.074±0.010 


0.137 


0.5+0.1 


0.113+0.014 


0.10 


1.1+0.1 


<0.67 






166.103 


14.8 


0.088±0.018 


0.204 


0.4+0.1 


0.114+0.021 


0.27 


0.4+0.1 


<0.10 






231410 


194 


<0.11 






0.325+0.038 


0.59 


0.5+0.4 


0.109+0.031 


0.182 


0.6+0.2 


275.671 


20.9 


<0.20 






<0.16 






<0.26 






276.060 


314 


<0.13 






0.092+0.051 


0.14 


0.6+0.4 


< 0.21 






CH2DOH 


89.408 




0.029±0.002 


0.053 


0.5+0.1 


0.040+0.007 


0.034 


1.1+0.2 


0.021+0.004 


0.057 


0.3+1.0 


CH3OD 


110.189 




< 0.04 






0.078+0.010 


0.073 


1.0+0.1 


0.025+0.013 


0.050 


0.5+0.3 



